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ABSTRACT 



As a part of a Euro-VO research initiative, we have un- 
dertaken a programme aimed at studying the scale length 
of 54909 Sa-Sd spiral galaxies from the SDSS DR6 cat- 
alogue. We have retrieved u, g, r, i, z-band images for 
all galaxies in order to derive the light profiles. We also 
calculate asymmetry parameters to select non-disturbed 
disks for which we will derive exponential disk scale 
lengths. As images in different bands probe different 
optical depths and stellar populations, it is likely that a 
derived scale length value should depend on waveband, 
and our goal is to use the scale length variations with 
band pass, inclination, galaxy type, redshift, and surface 
brightness, in order to better understand the nature of spi- 
ral galaxies. 
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1. INTRODUCTION 



The exponential scale length of a galaxy disk is one of 
the most fundamental parameters to determine its mor- 
phological structure, as well as to model its dynamics, 
and the fact that the light distributions are exponential 
makes it possible to contrain the formation mechanisms 
(Freeman 1970). The scale length determines how the 
stars are distributed throughout a disk, and can be used 
to derive its mass distribution, assuming a specific M/L 
ratio. Ultimately, this mass distribution is the primary 
constraint for determining the formation scenario (e.g., 
Dutton 2008, and references therein), which dictates the 
galaxy's evolution. As the galaxy evolves substructures 
such as bulges, pseudo-bulges, bars, rings, and spiral 
arms may build up, which in turn considerably change the 
morphology of the host disks (e.g., Combes & Elmegreen 
1993). Analytic disk formation scenarios (e.g., Lin & 
Pringle 1987) predict that in cases where angular momen- 
tum is conserved, the disk scale length is determined by 



the the angular momentum profile of the initial cloud, and 
the scale length in a viscous disk is set by the interplay 
between star formation and dynamical friction (e.g., Silk 
2001). These processes form the basis of a galaxy's grav- 
itational potential and the strength of gravitational pertur- 
bations, the location of resonances in the disk, the forma- 
tion and evolution of spiral arms and bars, and the dy- 
namical feeding of circumnuclear starbursts and nuclear 
activity (e.g., Elmegreen et al. 1996; Fathi et al. 2008). 

Photometrically, the scale length is derived by az- 
imuthally averaging profiles of the surface brightness 
which is in turn decomposed into a central bulge and an 
exponential disk, accounting for other components such 
as bars and rings. 

As images in different bands probe different optical 
depths and stellar populations, it is likely that a derived 
scale length value should depend on waveband. Dusty 
disks are more opaque and often deliver larger scale 
length values in bluer bands when compared with red 
and/or infrared images. Similar effects can also be caused 
by the stellar populations. These observational effects 
thus not only give us insights about the disks that we 
are studying, but also need to be quantified for a better 
comparison between different data sets and galaxy types. 
Both the effects of stellar populations and dist extinction 
have been subject to much discussion over the years (e.g., 
Simien & de Vaucouleurs 1983; Valentijn 1990; van Driel 
et al. 1995; Peletier et al. 1995; Beckman et al. 1996; 
Prieto et al 2001; Graham & de Blok 2001; Mac Arthur 
2003; Cunow 1998, 2001, 2004). A detailed and exten- 
sive analysis of the dust effects has also been presented 
for a few tens of galaxies in Holwerda (2005) and subse- 
quent papers, however, as noted by Peletier et al. (1994) 
and van Driel et al. (1995), the scale length alone in 
different band passes in small sample cannot be used to 
break the age/metallicity and dust effects. Investigating 
the scale length variation as a function of inclination for 
large numbers of galaxies is necessary to distinguish be- 
tween the dust and population effects. 
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Table 1. Sample sizes for the work mentioned in the text 
(in alphabetic order). 



Reference 



Number of galaxies 



Cunow (1998, 2001, 2004) 
Graham & de Blok (2001) 
MacArthur et al. (2003) 
Peletieret al. (1995) 
Prietoetal. (1996) 
Simien & de Vaucouleurs (1983) 
van Driel et al. (1995) 



14, 60, 39 
120 
121 
37 
15 
98 
55 



The common denominator in all the previous studies is 
the roughly comparable sample sizes. Most studies have 
so far analysed individual galaxies, or samples contain- 
ing a few tens of galaxies (see table 1). This is not to 
be mistaken with the number of great results from the 
Sloan Digital Sky Survey (SDSS) studies in the last years, 
but these works have not studied the astrophysical effects 
mentioned here. We have undertaken a programme that 
aims at quantifying how the disk scale length varies with 
band pass, inclination, galaxy type, redshift, and surface 
brightness. We have searched the entire SDSS Data Re- 
lease 6 (DR6) data set and have selected 54909 spiral 
galaxies suitable for our analysis. Here, we present a 
description of our study along with some preliminary re- 
sults. 



2. SAMPLE SELECTION 



The sample was selected by searching the SDSS DR6 
catalogue and cross matching with the LED A[[] catalogue 
(Paturel et al. 2003) to retrieve Hubble classifications. As 
the SDSS provides a number of morphological as well 
as kinematic parameters, which we use to constrain our 
sample against biases. Our first requirements ensure that 
for each galaxy we have reliable redshift measurement, 
low extinction, and number of pixels sufficient to derive 
a light profile with a good coverage of the disk region. 
Moreover, we decided to ensure that our galaxies do not 
contain edge-on or face-on galaxies to avoid selection ef- 
fect problems. We first retrieved tabular data for all SDSS 
galaxies for which excellent image quality is delivered, 
are larger than 30 pixels, spectroscopic redshifts are avail- 
able, have extinction Ay < 1.0 mag, have inclination 
15° < i < 70°. We retrieved a total of 475408 galax- 
ies, and first investigated the smallest numbers of pixels 
needed to resolve the disk. We found that a minimum of 
70 pixels are needed, thus removed all galaxies with ma- 
jor axis (in r-band) smaller than 70 pixels (28 arcsec) to 
ensure that the images cover the disk region. 

We made use of LEDA by first retrieving the entire cat- 
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Figure 1. Top: Right ascension and declination distribu- 
tion of the SDSS sample fulfilling the first sorting criteria 
(i.e., galaxies with spectroscopic redshift, larger than 10 
pixels, extinction smaller than Ay — 1.0, and with incli- 
nation between 15° and 70° ), Middle: The total LEDA 
catalogue for which the LEDA services provide a hubble 
classification number (0=S0a, l=Sa, 2=Sab, etc.), and 
Bottom: Our final sample fulfilling all our selection cri- 
teria described in the text. 



' http://leda.univ-lyon 1 .67 
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Figure 2. Distribution of some key parameters (magnitude, Galactic extinction, major axis, redshift, velocity dispersion, 
and morphological type) retrieved from the SDSS and LEDA database. 



alogue. As this service provides a numeric Hubble clas- 
sification parameter, we could easily select all the spiral 
galaxies, which we later cross-correlated with the SDSS 
catalogue. We found a total of 54909 Sa-Sd spiral galax- 
ies (see Fig. 1), for which all the morphological and spec- 
troscopic parameters from SDSS and LEDA were stored, 
and u, g, r, i, z-band images were to be downloaded. It 
should be noted that at this stage, we are unable to de- 
termine whether the galaxies in our sample are isolated 
or disturbed systems, as this information is not provided 
by any of the catalogues we have used. We make this 
distinction using the asymmetry parameter described in 
Conselice (2003). 

The first question that rises at this point is the fact that 
SDSS delivers the disk scale length as well as de Vau- 
coulers effective radius for each galaxy (in all bands), and 
that these values could be used to carry out or analysis. 
In Fig. 2, we show that the values provided by the SDSS 
team show anomalies that are beyond our satisfaction for 
carrying out our analysis. The plot shows a strange "clus- 
tering" of the effective radii and scale lengths around 
some numbers, the source for which we cannot find. We 
thus decide to re-calculate the scale lengths. 

Various Virtual Observatory (VO) methods were investi- 
gated to perform the download of the SDSS images. The 
Sky View was chosen for this task. This service has the 
advantage of being able to create fits files centred at a 
given sky coordinate and with a pre-specified size. The 
image size is an important parameter for achieving a re- 
liable an accurate sky subtraction, thus we require that 




SDSS disk scale length 



Figure 3. De Vaucouleurs effective radius (y-axis) versus 
disk scale length (x-axis) provided by the SDSS team for 
all 54909 galaxies in our final sample. Although these 
numbers are good first guesses, the odd clustering of the 
data points lead us the to conclusion that we re-derive the 
scale lengths. The insets illustrate the strange clustering 
of the points ( coloured data points) provided by the SDSS. 



2 http://skyview.gsfc. nasa.gov/ 
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Figure 4. Illustrating the fits to the different regions of the light profile. Fitting different regions results in different scale 
length values, and as shown in the bottom right panel, fitting over the entire galaxy reproduces the value provided by 
the SDSS, which is illustrated by the "blue" line (arbitrarily shifted (in y-axis). Printed at the top of each panel, is the 
corresponding scale length value. 



the images are 900 x 900 pixels for the sky region to be 
sampled for all galaxies. Moreover, SkyView is able to 
re-scale the image backgrounds to the same level, hence 
correcting for differences between the SDSS plates. We 
have experienced that a linear image download could be 
a tedious process. With an image size of 3.2 MB/image, 
and at a constant rate of 0.2 MB/s, the download would 
require at the very best and with a continuous connection 
f» 51 days. We have therefore carried out the download 
using multiple parallel data requests to the SkyView ser- 
vice using the python scripting language, and utilising 3 
individual computers. This allowed us to download the 
full sample over approximately 13 days. 



3. DERIVING THE LIGHT PROFILE 

We derive the disk scale length using standard IDL rou- 
tines and make use some important parameters included 
in the SDSS information in order to constrain galaxy ge- 
ometry as well as the location of the sky region. These 
are isop, isoa, isophi, and for consistency, we use these 
quantities in r-band. The procedure carries out the fol- 
lowing steps: 

• Reading the image, and calculating the asymmetry 
parameter A = J^(7 — Ixso) I X} I> where I is the 
galaxy image and Iiso is the image rotated by 180 
degrees around the galaxy centre. 



• Selecting the sky region as the ellipse encompass- 
ing the region between 2.0 * isoa and 2.5 * isoa- 
The mean value of this region, using Tukey's bi- 
weight mean formalism described in Mosteller & 
Tukey (1977), is used for the sky subtraction as well 
as setting the background level. 

• To remove nuisance stars and point sources from 
the image, we extract point sources with SExtractor 
(Bertin & Arnauts 1996), but storing pixels belong- 
ing to all point sources that are larger than 4 pixels 
and more than 3 a above the background level. All 
pixels belonging to these sources are then masked 
out. 

• Using the isop, isoA, isophi parameters from 
SDSS, we then section each galaxy into 2-pixels 
wide ellipses oriented at isophi and with minor-to- 
major axis ratio b/a = isop /isoa- The mean sur- 
face brightness value within each ellipse is calcu- 
lated, to compile the galactocentric light profile for 
each galaxy image. 

• As experience has shown that light profile decompo- 
sition in a disk together with any other component 
(e.g., Fathi et al. 2003) introduces complications 
that are not necessary for the nature of our analy- 
sis. We thus derive the disk scale length simply by 
fitting an exponential profile to the disk region of 
each galaxy. We determine the region of interest by 
empirically fitting eight different regions covering 
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the range between 15%-30% and 85%-115% of the 
isoA parameter. This procedure means that we are 
simply cutting out the central regions of the galaxies 
where bulges and strong bars are expected. Fig. 4 
shows the result of such a test where we find that 
(in many cases) we are able to derive a scale length 
comparable to the value from the SDSS catalogue. It 
should be noted that this plot is only for one galaxy, 
and for regions outside what is noted here, and we 
have tested more regions than the eight regions de- 
tailed above here. 



Assuming that the r and z-band images probe similar stel- 
lar populations and dust content, we are able to use the 
correlation of the scale length values between these bands 
to find the optimal region for deriving the scale lengths. 
Fig. 5, illustrates for a preliminary and randomly selected 
sub-sample of 1315 galaxies, that using a small outer re- 
gion is not the best way to derive the scale length, but for 
this sub-sample is seems that fitting the galaxy light pro- 
file over 1.2 * isoA gives the best result. We plan to carry 
out this test for the entire sample. 



4. CURRENT STATUS 



We have currently downloaded the u, g, r, i, z images for 
all 54909 galaxies, and we are in the process of calcu- 
lating the scale length using different regions from each 
image. With the data at hand, we will be able to first 
remove highly asymmetric galaxies to minimise the use 
of disturbed disks. We will then derive the disk scale 
lengths for all the "isolated" disk galaxies, and statisti- 
cally explore how this parameter changes as a function 
of inclination, band, redshift, etc. Moreover, we plan to 
cross-correlate our sample with the Two Micron All Sky 
Survey to further explore the scale lengths also in J, H, K 
bands and to further explore the dust and stellar popula- 
tion effects. 



Figure 5. Scale lengths for a sub sample of 1315 ran- 
domly selected galaxies in u, r, i, and z-band. Here, fits 
to three different disk regions are plotted to show that fit- 
ting some outer region is not necessarily the best strategy. 
The relatively large scatter in the top and bottom panels 
also demonstrate that the images are not very deep in u 
and z bands. At the bottom right corner of panel, the 
fitted regions are printed, where 0-120% means that the 
light profile between galactocentric radius ofr = 0*isoA 
and r — 1.2 * isoa was fitted by one single exponential 
profile, etc. . 
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